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ABSTRACT 

: 

m . 

. Gravitational waves emitted by chirping supermassive black hole binaries 

could in principle be used to obtain very accurate distance determinations. Pro- 
vided they have an electromagnetic counterpart from which the redshift can be 
determined, these standard sirens could be used to build a high redshift Hubble 
diagram. Errors in the distance measurements will most likely be dominated by 
gravitational lensing. We show that the (de)magnification due to inhomogeneous 
+3 . foreground matter will increase the scatter in the measured distances by a factor 

d ■ ~ 10. We propose to use optical and IR data of the foreground galaxies to mini- 

£> . mize the degradation from weak lensing. We find that the net effect of correcting 

^ . the estimated distances for lensing is comparable to increasing the sample size 

by a factor of three when using the data to constrain cosmological parameters. 

Subject headings: cosmological parameters — gravitational lensing — gravita- 
tional waves 



1. INTRODUCTION 

The distance-redshift relation for electromagnetic standard candles, such as Type la 
supernovae, has been crucial to our understanding of the energy content of the Universe. 
With the advent of laser interferometric gravitational wave detectors, such as LIGO0 and the 
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planned space mission lisa| new windows to the Universe will be opened. Binary systems 
with shrinking orbits, emitting gravitational waves just before coalescence, are expecte d to be 
standard si rens, the gravitational wave analogy of electromagnetic standard candles (jSchutz 
19861 . 120021 ). Since the frequency of the gravitatio nal waves ris es when the orbit changes, 
these systems are also known as chirping binaries (ISchutzlll986l ). Laser interferometers can 
provide crude estimates of the position of the gravitational wave source in the sky. If a 
source has an observable optical counterpart, the posi tion of the si ren could be pinpointed 
and the distance measurement significantly improved (jSchutzl 119861 ). 



Accurate redshifts of electromagnetic standard candles are usually easier to obtain than 
accurate distance measurements. For gravitational wave standard sirens the situation is 
reversed: accurate distances, but no redshifts, can be obtained. An electromagnetic counter- 
part could thus not on ly increase the accuracy in the distance measurable by LISA to < 1% 
(IHolz fc Hughesll2005l ) but also provides the redshift. Detection of gravitational wave sources 
with electromagnetic counterparts th us have the p otential to probe the distance-redshift 
relation with unprecedented precision (ISchutzlll986l ). For a recent discuss ion on possible ob - 
servable electromagnetic counterparts of massive black hole binaries, see iDotti et al.l (120061 ) 
and references therein. 

The coalescence of supermassive black holes is associated with mergers of proto-galaxies 
or galaxies, many of which are expected to occur at high redshift. Chirping supermassive bi- 
nary black holes thus promise to be probes of the very distant Universe. The exact redshift 
distribution of supermassive binary black holes is, however, poorly known and the esti- 
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20061 ). While cosmological parameters can be derived from a Hubble diagram of "standard 
sirens" at high z, severe parameter degeneracies result unless there are also data points at 
low redshifts. Since the distances which could be obtained from chirping supermassive black 
hole binaries are absolute, they could be combined with abs olute distances m easured using 
other low redshift probes, e.g. chirping neutron star binaries ( IDalal et al.ll2006l ). At high red- 
shifts, distance mea surements will be affected by gravitational deflection of the radiation by 
foreground matter (lMarkovidll993l ). An additional scatter of a factor ~ 10 is thus expected 
due to gravitational lensing (see Section @J. 

In principle, due to flux conservation, the effects of gravitational lensing are expected 
to average out if a large number of standard candles or sirens are observed. However, 
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the large uncertainty in the rate, coupled with the fact that it might be difficult to find 
and observe electromagnetic counterparts, suggests that a scenario where only a handful of 
standard sirens with redshifts can be measured by LISA is not unlikely. The scatter due to 
gravitational lensing will then limit the usefulness of the data severely. 

In this pa per, we propose a method to compensate for the degradation due to lensing, at 



least partially. iGunnarsson et al.l (120061 ) showed that properties of foreground galaxies can be 
used to estimate the magnification accurately enough for corrections of individual sources to 
be useful. The error box in the sky provided by LISA, will probably be extensively observed 
in many different wavelength bands in the search for an electromagnetic counterpart. This 
data could be used to estimate the magnification of the siren due to gravitational lensing 
by foreground matter. For reasonable detection rates of chirping supermassive binary black 
holes, we show that correcting for gravitational lensing will be comparable to increasing the 
sample size by a factor three. The corrections will thus be crucial when using the data to 
constrain cosmo logical parameters such as the Hubble constant, the matter density and dark 
energy properties. The cosmology used throughout this paper is the concordance flat ACDM 
model, [h,Qu,wo] = [0.7,0.3,-1]. 



2. CORRECTING GRAVITATIONALLY LENSED SOURCES 

In this section we describe how standard sirens can be corrected for gravitational lensing. 



2.1. Computing Corrections 

As noted, gravitational waves emitted from high redshift standard sirens might be grav- 
itationally deflected by matter along the line of sight. The magnification factor, /i, describes 
the amplification (ji > 1) or de-amplification (p < 1) of a source due to gravitational lensing 
as compared to a universe with homogeneously distributed matter. If is the luminosity 
distance in a homogeneous universe (which is what is used when estimating cosmological 
parameters), the observed luminosity distance, -D£ bs , affected by gravitational lensing, is 



Gunnarsson et al.l (120061 ) showed that observed properties (e.g. position, spectroscopic 
or photometric redshift, and luminosity) of intervening galaxies along the line-of-sight of 
the source can be used to estimate the magnification factor with sufficient accuracy for 
corrections to be useful. The magnification factor can be computed using halo profiles 
depending on the mass of the galaxy halo, which can be obtained from the galaxy luminosity 
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using Tully-Fisher and Faber- Jackson relati ons. This method ha s been successfully applied 
to, e.g. 33 supernovae in the GOODS fields (IJonsson et al.ll2006l ). The luminosity distance 
corrected for gravitational lensing using the estimated magnification factor, // est , is 

or = ^Df s = (fi/^ s r l/2 D L . (i) 

Whether the corrected distance is closer to the obtained in a homogeneous universe 
than the observed distance depends on how close the quantity /i//i es t is to unity, i.e. how 
well we can estimate the magnifications. This in turn depends on the completeness and 
accuracy of the foregrou nd galaxy information. We have used the method presented in 



Gunnarsson et al.l (120061 ) to obtain the precision to which /z est can be estimated for sources 



out to z — 4. The distribution of /////est at a specific redshift was found by computing \i 
and //est for a large number of simulated lines-of-sight. Galaxies along each line-of-sight were 
simulated and the magnification computed using a fiducial model. This value was taken as 
the true magnification /t. Since the mean of a simulated probability distribution function 
(PDF) of // is not exactly unity, each PDF of // was normalized to ensure a mean of unity. 
A value of the magnification representing /t est was also computed assuming a limiting galaxy 
magnitude and a reasonable error budget including the effects of using the wrong lens model 
(for 50% of the lenses a SIS-profile was used instead of a NFW-profile) and uncertainties 
in redshifts and positions as well as intrinsic scatter and accuracy in the Tully-Fisher and 
Faber- Jackson relations. Since new powerful instruments for electromagnetic observations 
should be available when LISA can be expected to be operational, our simulations assume 
galaxy optical limiting magnitude of / = 28. In these simulations we have only considered 
the gravitational deflection due to galaxies and neglected any effects of large scale structures. 



2.2. Distance Probability Distribution Functions 

The accuracy to which the luminosity distance to chirping supermassive binary black 
ho les could be measured by LISA, neglecting the effects from lensing, has been investigated 



by lHolz fc Hughes! (120051 ). In their most optimistic scenario, where an electromagnetic coun- 
terpart was assumed, the fractional error in the distance measured for sirens at redshift z = 1 
and z = 3 was found to be 5D^/D^ ~ 0.1% and SD^/D^ ~ 0.5%, respectively. We use these 
numbers and linear interpolation and extrapolation to obtain the fractional errors used in 
subsequent calculations. The PDF describing the dispersion in the distance obtained from 
an unlensed standard siren is assumed to be a Gaussian. 

Supermassive binary black holes thus have the potential to be almost perfect (absolute) 
standard sirens and the gain from being able to minimize any additional scatter, e.g. lensing, 
can be very large. 
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The PDF of the observed luminosity distance D£ bs = /i~ 1//2 DL is given by the convolution 
of the PDFs of fi^ 1 ^ 2 and D^. Both PDFs depend on redshift and, to some extent, on 
cosmological parameters. Corrected luminosity distances are expected to be distributed 
according to the PDF given by the convolution of the PDFs of (/i//i est ) _1//2 and D^. The 
cosmology dependence on the magnification factor is weak and is therefore neglected. 

Figure [1] shows the PDFs of luminosity distances measured from unlensed, lensed, and 
corrected standard sirens at different redshifts. Solid lines indicate the luminosity distance 
PDF for unlensed sirens. Gravitational lensing introduces asymmetry and additional disper- 
sion to the observed luminosity distance PDF as can be seen from the dotted curves. If the 
observed distances are corrected for gravitational lensing, both dispersion and asymmetry 
can be significantly reduced. Corrected PDFs are indicated by dashed lines. 

Although the mean magnification factor is unity, the mean of /z -1 / 2 might be different 
from unity depending on the moments of the PDF of /i. The mean luminosity distance of a 
lensed or corrected standard siren might thus be biased compared to an unlensed siren. This 
bias can be circumvented if D^ 2 is considered instead of D L (see Section Hj). For a lensed 
standard siren at z = 1 the mean luminosity distance is overestimated by ~ 0.1% and this 
effect increases with redshift to ~ 0.6% at z = 4. Correcting for gravitational lensing reduces 
the bias at z < 3, but increases the bias at z = 4 where the mean of the corrected luminosity 
distance is overestimated by ~ 1%. However, since we expect only one or a few standard 
sirens at each redshift, the location of the mode is more important than the mean. As can 
be seen in Figure [Tj, the mode of the corrected luminosity distance PDF is always closer to 
the true unlensed distance than the mode of the observed (lensed) luminosity distance. 

After correcting for gravitational lensing, we are left with an absolute standard siren 
with an almost Gaussian dispersion of SD L /D L ~ 1% and 5D L /D L ~ 3% at redshift z = 1 
and z = 4, respectively. 

2.3. Distance Error from Position Uncertainty 

The accuracy to which the position of the standard siren could be determined affects 
not only the distance determination, but also the uncertainty of the estimation of the mag- 
nification factor. If the position of the siren could be narrowed down to the host galaxy 
only, the dispersion in the corrected luminosity distance could be significantly increased. An 
uncertainty of 3 arcseconds in the position of the source increases the dispersion in the cor- 
rected luminosity distance with 10% and 60% at redshift z = 1 and z — 4, respectively. In 
the following, we assume that the position of the source is pinpointed with enough precision 
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for this source of error not to dominate. Note, however, that this corresponds to arcsecond 
resolution at z = 4 which could be an optimistic assumption. 



2.4. Groups, Clusters, and Dark Halos 



In the previous sections, we have assumed a one-to-one relation between luminous galax- 
ies and dark matter halos. If this assumption is relaxed, our ability to predict lensing 
magnifications could potentially be severely impaired. This could be caused by, e.g. com- 
pletely dark halos whitout any galaxies, and dark matter halos hosting groups or clusters 
of galaxies. We have investigated the size of this effect by simulating lines-of-sight using 
public ly available halo and ga laxy catalogs constructed from the results of N-body simula- 
tions (IKauffmann et al.lll999l ). Particles in the simulations which were not bound to any 



halo w ere treated as if they belonged to a smooth matter component (see iGunnarsson et al. 



(120061 ) for details). To compute the true magnification factors, data from the halo cata- 
logs, including halo masses, were used. When the magnification factor was estimated, data 
from the galaxy catalogs wer e used and halo m a sses w ere obtained from galaxy luminosities 
following the pr escription in IGunnarsson et al.l (120061 ). However, we used a different Tully- 



Fisher r elation (jGiovanel 



used by IKauffmann et al 



i et al 



19971 ) in order to match the normalization of the models 



(119991 ). In the calculations of the magnification factors all halos 



were parametrized as truncated NFW-profiles. 

We have studied the effects of groups, clusters, and dark halos on our ability to predict 
the magnification of a source at z — 1.5. For a source at this redshift the dispersion, in 
terms of the root mean square deviation, of the PDF of \x is ~ 7% and the corresponding 
dispersion in the luminosity distance PDF is ~ 3%. If the magnification could be recovered 
exactly, the dispersion in /x//x es t would vanish. 



The effect of dark halos is partly taken into consideration in IGunnarsson et al.l (120061 ) 
where galaxies fainter than the magnitude limit play the role of completely dark halos. 
However, to make an extremely conservative estimate, we use simulated lines-of-sight from 
N-body results where o nly the most massive, but not all, of the halos contain galaxies 
( IKauffmann et al.lll999l ). The contribution to the dispersion in /i//i e st due to the ~ 90% of 
the halos in the simulated lines-of-sight which are not associated to any galaxy and thus 
cannot be modelled, is ~ 3%. 

Approximately 2% of the dark matter halos host more than one galaxy, corresponding to 
groups or clusters of galaxies. Even though this number is low, since these halos are also the 
most massive ones, the effect of poor modelling of groups of galaxies could potentially add a 
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sizeable dispersion. Our simulations show that these halos add a dispersion of ~ 3% to the 
dispersion of /x//i est . However, since groups and clusters of galaxies could be identified, better 
modelling of the halos of these systems might be used, which would reduce the dispersion. 
For lines-of-sights in the outskirts of clusters even lensing maps constructed from multiple 
image systems might be used to estimate the magnification. 

If a ~ 50% scatter in the relation between galaxy luminosity and mass, completely 
dark halos, and poor modelling of galaxy groups are taken into account, the dispersion 
in the magnification factor can be reduced from ~ 7% to ~ 4%. This translates into an 
improvement in the dispersion in the luminosity distance from ~ 3% to ~ 2%. Such a 
reduction of the dispersion corresponds to increasing the sample size by almost a factor 
3. However, the mean of the corrected luminosity distance is underestimated by ~ 1%, 
considerably worse than the bias of ~ 0.2% toward higher mean in the case of observed 
(lensed) luminosity distance. 

A biased relation between galaxy luminosity and halo mass could severely impair our 
capability to correct for gravitational lensing. If all masses were underestimated by ~ 50%, 
the dispersion in the magnification factor and luminosity distance could be reduced to ~ 
5% and ~ 2%, respectively. Since the the mean of the corrected luminosity distance is 
underestimated by only ~ 0.5% in this case, such a bias in the luminosity-mass relation 
would not significantly hamper or ability to correct for lensing. 

If, on the other hand, all masses were overestimated by ~ 50%, correcting for gravita- 
tional lensing would instead increase the dispersion in /i//i es t to ~ 10%, ruining our attempt 
to correct for gravitational lensing. Since we require the total mass in halos to be < f^M, a 
too large bias toward higher halo masses would be discovered. Furthermore, any substantial 
bias in the luminosity-mass relation should give worse residuals in the Hubble diagram after 
correcting data for lensing than before. 



In the so called realistic scenario outlined in iGunnarsson et al.l (120061 ). where several 
effects are taken into account, the dispersion in /i//i es t for a source at z — 1.5 is ~ 3%. The 
discussion above lead us to the conclusion that in order for lensing corrections to be useful, 
the dispersion in /i//i cst must be < 5%, leaving us with an error budget of ~ 4% for dark 
halos and groups and clusters of galaxies. Since our extremely conservative analysis yields 
< 4%, we are confident that this number is within reach. 
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3. IMPROVING COSMOLOGICAL PARAMETER ESTIMATION 

To illustrate the improvement if corrections for gravitational lensing are applied, we 
have fitted cosmological parameters to simulated data sets and plotted the corresponding 
confidence contours. 



3.1. Simulated Data Sets 



Our fiducial supermassive binary black hole data set consists of 10 luminosity distance 
measurements to standard sirens evenly distributed in the redshift range 1 < z < 4 with 
errors as described in Section [2T21 (high redshift sample). In addition, chirping neutron star 
binaries might provide gray itational wave standard sirens with optical counterparts at low 
redshifts ( IDalal et al.ll2006l ). We assume 5 additional distance measurements at low red- 
shifts, 0.04 < z < 0.12 (low redshift sample). The fractional err or in these distanc es were 
assumed to be given by SD^/D^ = Z?l/(1700 Gpc) adopted from IDalal et al.l (120061 ) . Since 
the improvement from correcting for lensing corresponds to a factor of three in sample size 
almost independently of the exact rate and redshift distribution of standard sirens (see Sec- 
tion H]), the fractional improvement in the confidence contours can be appreciated although 
the absolute number of well observed sirens is highly uncertain. 



3.2. Maximum Likelihood Technique 

Since the PDFs of lensed distances are expected to be asymmetric (see Figure [TJ, a 
maximum likelihood technique is used to fit cosmological parameters to the simulated data. 
Given a data set of redshifts and luminosity distances, {zi, Di}, obtained from gravitational 
wave standard sirens and their electromagnetic counterparts, cosmological parameters can 
be estimated by maximizing the logarithm of the likelihood, 

log£(0) = £log/(A;*,0), (2) 

i 

where 9 represents the cosmological parameters whose values we want to estimate. We fit 
parameters to data sets consisting of unlensed (D = Dj), lensed (D = /i -1 / 2 ^), and cor- 
rected (D = (/i//i est ) _1 / 2 .D L ) luminosity distances. For each data set the corresponding PDF, 
f(D;z,9), is used. Since the effects of gravitational lensing are redshift dependent, PDFs 
for lensed and unlensed luminosity distances are redshift dependent too. Low redshift sirens, 
with negligible magnification due to lensing, are assumed to be described by a Gaussian 
PDF. 
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In the cosmology fits, a model for luminosity distances in a flat universe is used, 

f z dz' 

D h (z;h,Qu,w ,w a ) = (1 + z) . (3) 

Jo u \ z ) 

This model depends on the Hubble constant, Hq = 100/i km s™ 1 Mpc -1 , the m atter density , 



Qu, and the equation of state parameter of dark energy parametrized by (e.g. iLinderl 120031 ) 
w(z) = w + w & zj (1 + 2;), via the Hubble parameter, 

H 2 (z) = Hi [Q M (1 + zf + (1 - fi M )(l + z )3(i+-o+»a) e -3^ 2 /(i+,)j _ 

Note that chirping binaries, in contrast to e.g. Type la supernovae, are absolute distance 
indicators, i.e. we don't have to marginalize over any non-cosmological parameters. We have 
considered the effects on confidence level contours for both constant and time varying dark 
energy equation of state parameters. In the cases where the equation of state parameter was 
assumed to be constant, the parameter w a was put to zero and one of the three remaining 
parameters was either marginalized over or assumed to be exactly known. When a time 
varying equation of state parameter was considered, prior information on Qm and h was 
used. 



3.3. Effects on Cosmology Fits 

Figure [2] shows the results of the cosmology fits where a constant equation of state 
parameter was assumed, i.e. w a = 0. Each panel in the figure shows 68.3% confidence 
level contours (log L(9) = logL max — 2.30/2) for unlensed (solid contours), lensed (dotted 
contours), and corrected (dashed contours) data sets. Each column of panels in the figure 
shows the results of cosmological parameters fitted to different data sets. The left column 
shows results obtained for the high redshift sample. The middle column shows the effect of 
adding the low redshift sample. The right column of panels shows the results of cosmology 
fits where the third parameter was assumed to be exactly known. This gives an indication of 
the dependence of the size of the contours on the marginalization over and prior knowledge 
of the third parameter. 

Panel a shows confidence level contours in the [Qm> w o] -plane obtained for a simulated 
data set consisting of 10 high redshift standard sirens. Gravitational lensing degrades the 
data set and the cosmological results that can be obtained from it. However, the size of the 
confidence level contours can be significantly reduced if the data are corrected for lensing 
magnification. Distance measurements at low redshifts are much more useful in constraining 
h than distance measurements at high redshifts. Adding low redshift distance measurements 
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to the high redshift data set thus have effects similar to applying prior information on h. 
Panel b shows the improvement on the confidence level contours due to 5 additional distance 
measurements at low redshifts. Panel c shows the very tight constraints in the [Om > ifo]-pl ane 
which could be obtained if h was exactly known. 

The middle row of panels shows results similar to those presented in the top row of 
panels, but in the [h, f2 M ]-plane. Panel d shows the contours obtainable from 10 high redshift 
standard sirens. The improvement due to 5 additional low redshift standard sirens can be 
seen in panel e. A data set consisting of both high and low redshift distances could potentially 
give constraints on flu and h on the 10% level. Panel f shows the results for the case where 
Wo is assumed to be exactly known. Note that assuming that the da rk energy is accurately 



descr ibed by a cosmological constant, as current data favors (e.g. lHannestad fc Mortsell 



2004 ). standard sirens have the potential to give very strong constraints on the Hubble 



parameter, h, if we correct for gravitational lensing. 

The bottom row of panels show results in the [h, u>o]-plane. Panel g shows confidence 
level contours obtained from 10 high redshift sirens. The improvement due to 5 additional 
low redshift standard sirens can be seen in panel h. Panel i shows the results for the case 
where Qm is assumed to be exactly known. 

Results obtained from the high redshift data sets of lensed or corrected sirens suffer from 
severe degeneracies. The (one-sided) limits on Qm, wo and h which could be obtained from 
the high redshift sample alone, without any priors or low redshift distance measurements, 
would not be competitive compared to other measurements of these parameters. Priors or 
low redshift data seem to be necessary for a small set of high redshift standard sirens to be 
able to give competitive constraints on Qm, wo and h. 

Figure [3] shows the results in the [w , w a ]-plane. Solid, dotted, and dashed lines indicate 
the 68.3% confidence level contours obtained for unlensed, lensed, and corrected data sets. 
Gaussian priors, assuming h and I^m to be known with 3% accuracy, were used in the fitting 
procedure. Panel a shows the contours obtained from the high redshift data set. The effect 
of adding the low redshift data set is small and is shown in panel b. Since strong priors have 
already been applied, there is no dramatic improvement due to broken degeneracies. Panel c 
shows the contours obtained with both h and flu exactly known. The set of 10 high redshift 
sirens is evidently not sufficient to yield very strong constraints on wq and w a even in the 
case where the Hubble constant and the matter density are exactly known. Larger data sets 
than the ones considered here are thus needed for standard sirens to be able to provide very 
strong constraints on the evolution of dark energy. 



In Table [T] and [2] the numerical results are summarized. The 68.3% confidence level 
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intervals in the tables were obtained by finding the extremal values of the cosmological 
parameter 6 (where 9 is h, f2 M , w , or u> a ) fulfilling the condition \og L{6) = log L max — 2.30/2. 



4. FLUX AVERAGING 

For a large number of sources, all at the same redshift, the mean magnification is 
expected to be unity due to flux conservation. Since gravitational wave flux is inversely 
proportional to the luminosity distance squared, the mean of -DL 2 obs = I^D^ 2 equals D^ 2 . 
The mean of the observed luminosity distance is, on the other hand, expected to be biased 
by 0.1% — 0.6% as noted in Section 12.21 This bias is similar to the one which arises if 
magnitudes are used instead of fluxes. However, whether this bias is important or can be 
ignored depends on the sample size and the errors in each distance measurement. In the 
scenario studied in this paper, where the standard sirens are sparsely distributed in redshift, 
the position of the mode is more important than the mean and we can use .D£ bs instead of 
.Dl obs as our observable. 

Degradation by gravitational lensing can thus be cured by studying the mean of the lumi- 
nosity distance (or inverse of the luminosity distance squared if it is importan t to avoid a bi- 



ased mean) with uncertainty crg L , if the number of sources is large enough (IHolz fc Linder 



20051 ). We now estimate how large this number of supermassive binary black holes is. 



Figure H] shows the dispersion of the luminosity distance PDFs as a function of redshift, 
quantified in terms of the root mean square of the deviation from the mean, <Jd l - The 
increase in the dispersion, due to lensing, depends on the redshift. At z = 1 and z = 4 
the dispersion is increased a factor ~ 20 and ~ 10, respectively. If the distances obtained 
from the sirens are corrected, the dispersion can be reduced by a factor ~ 2. If a number 
of lensed sirens, all at the same redshift, would be observed, how many would be needed 
to compensate for the dispersion due to lensing? The dispersion of the mean value of the 
distances inferred from the standard sirens should be inversely proportional to the square 
root of the number of sirens, 

„■ flensed at— 1/2 /r\ 

a S L = v Dl n 1 (5) 
To restore the dispersion to the value corresponding to one unlensed siren, i.e. = 
^uniensed^ needed number of observed sirens would thus be, 

a t /lensed / unlensed \ 2 (a\ 

Since the dispersion of an unlensed siren increases much faster as a function of redshift than 
the dispersion of a lensed or corrected siren, the number of sirens needed to compensate for 
lensing is much larger at z — 1 (N ~ 500) than z = 4 (N ~ 70). Only in very optimistic 
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scenarios would we be able to observe such a large number of sources including optical 
counterparts. 

The number of lensed sirens needed to achieve a dispersion comparable to the dispersion 
of one corrected siren is 

at /lensed / corrected \ 2 o fr*\ 

N = (^Dl /^Dl ) ~ 3, (7) 

almost independent of redshift. Thus, unless the number of well-observed sirens is very large, 
by correcting for the effects of gravitational lensing, we are able to improve our data by an 
amount that is comparable to increasing the sample size by a factor of three. 



5. SUMMARY AND DISCUSSION 

Gravitational waves emitted by chirping supermassive black hole binaries measured by 
LISA could in principle provide very accurate distance determinations. If electromagnetic 
counterparts of these events could be found, from which redshifts could be determined, these 
standard sirens could be used to build a high redshift Hubble diagram. In principle, chirping 
binaries are almost perfect distance indicators, only limited by the precision of the position 
measurement. Thus, the need of controlling systematic errors on the distance measure is 
crucial to be able to use the standard sirens to their full potential. These errors will most 
likely be completely dominated by gravitational lensing. 

Since the number of expected events is uncertain and whether optical counterparts 
exist or can be observed is even more unclear, we have reasons to believe that the number of 
events which could be put in a Hubble diagram is low. Thus, we can not be certain that the 
degradation due to gravitational lensing can be overcome by large statistics. From Figure [2] 
and (31 it is obvious that correcting for gravitational lensing can be crucial for the competitive 
use of gravitational wave standard sirens in cosmology. Note also that the actual number 
of observed sirens with electro magnetic counterparts at very high redshift (z > 4) might be 



substantial (IDotti et al.ll2006l ). in which case the effects from gravitational lensing will be 
even larger. 

Low redshift points in the Hubble diagram (where the luminosity distance mainly is 
sensitive to the Hubble constant) are important to break some of the degeneracies on the 
cosmological parameters occurring for high redshift luminosity distances, but no supermas- 
sive black hole binaries are expected at low redshift. The situation is thus very similar to 
the case of using gamma ray bursts as standard candles where, apart from the difficulty of 
calibrating the candles, th e lack of low redshift sources weakens their power to constrain cos- 



mological parameters, see iMortsell fc Sollermanl (120051 ). Since chirping supermassive black 
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hole binaries are absolute standard candles however, their usefulness for cosmology could 
be increased if they are combined with other absolute distance determinations in the local 
universe. 

Besides constraining the distance-redshift relation, standard sirens have other potential 
applications in cosmology. Since the dispersion from gravitational lensing depends on the 
matter distribution in the universe, it could in principle be used to constrain cosmological 
parameters. Equivalent to the use of galaxy shear measurements, the lensing magnifica- 
tion PDF can provide information on dg and f^M- The feasibility of the method has been 



studied in the context o f futu re Type la supernova observations in ICooray et al.l (120061 ) and 



Dodelson &: Vallinottol (120051 ). In this paper, we only note that standard sirens have the 



potential to be even better probes of the cosmic magni fication, due to the much smaller (at 



least a factor of ten) intrinsic scatter in the luminosity (ICooray et al.ll2006l ). Also, the large 
distances to the standard siren sources makes it possible to study the matter distribution in 
the universe at high redshifts. However, whether this potential can be realized depends upon 
the very uncertain detectability rate and the possibility to secure redshifts for the sources. 

Even if the number of observed sirens turns out to be low, the small intrinsic dispersion 
combined with the fact that chirping binaries are absolute standard candles allows us to 
measure the magnification for individual sources and thus probe the matter content, i.e., the 
convergence k along individual lines-of-sight. This could yield important information on the 
matter distribution in clusters and galaxies. Specifically, it could be used to break the mass 
sheet degeneracy that plagues galaxy shear investigations. 

In conclusion, we note that although detection rates are very uncertain, LISA could 
in principle provide high precision distance measurements to very distant objects in the 
Universe, the error of which will be dominated by gravitational lensing. Correcting for these 
errors is comparable to increasing the sample size by a factor of three which decreases the 
error contours for cosmological parameters substantially. Depending on the exact redshift 
distribution of the sirens, this could yield competitive constraints on, in particular, the 
Hubble constant and the matter density, but perhaps also on dark energy properties. 
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Fig. 1. — Probability distribution functions (PDF) for luminosity distances obtainable from 
unlensed, lensed, and corrected standard sirens at different redshifts. The fractional uncer- 
tainty to which the luminosity distance to an unlensed siren could be measured by LISA is 
indicated in each panel. 



-17- 



o 
- 1 

-2 
-5 
-4 



0.0 0.1 



0.2 0.3 



0.4 0.5 



a 



O.b 
0.4 

0.3 

0.2 

0. 1 
0.0 
0.5 



0.6 



0.7 0. 
h 



0.9 1.0 




- I 

-2 

-5 

-4 



0.5 0.6 



0.7 0. 
h 



0.9 1.0 



0.0 ~~ 
-0.5 - 

5° -1.0 " 

-2.0 _ 



0.3 ~ 
0.4 - 

a 0.5 - 
0.2 - 



u. 



0.0 ~ 
-0.5 - 

s° -1.0 - 
-2.0 _ 



0/ 0.2 



0.5 0.6 



0.5 0.6 



0.3 



0.7 
h 



0.7 

h 



0.4 0.5 



0.8 0.9 



0.8 0.9 



-0.6 ~ 
-0.8 - 

-1 .0 - 
-1 .4 _ 



0.20 0.25 0.30 0.35 0.40 

0.40 I ' ' ' 7~ 
0.35- \\ 

0.50 - "% ; , 

0.25 - 

0.20 I , , , 

0.60 0.65 0.70 0.75 0.80 
h 



-0.6 ~ 
-0.8 - 

-1 .0 - 
-1 .4 _ 



0.60 0.65 0.70 0.75 0.80 
h 



Fig. 2. — Confidence level contours at the 68.3% confidence level for different cosmological 
parameters and different simulated data sets. Solid, dotted, and dashed contours corresponds 
to unlensed, lensed, and corrected data sets, respectively. Panel a: High redshift data set. 
Panel b: High and low redshift data sets. Panel c: High redshift data set but h is known 
exactly. Panel d: High redshift data set. Panel e: High and low redshift data sets. Panel f: 
High redshift data set but wq is known exactly. Panel g: High redshift data set. Panel h: 
High and low redshift data sets. Panel i: High redshift data set but Qu is known exactly. 
Note the different scales in the figures. 
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Fig. 3. — Confidence level contours at the 68.3% confidence level in the [wq, w a ]-plane for 
different data sets. Solid, dotted, and dashed contours corresponds to unlensed, lensed, and 
corrected data sets, respectively. Panel a: High redshift data set. Panel b: High and low 
redshift data sets. Panel c: High redshift data set but h and f2 M are known exactly. Gaussian 
priors (3% accuracy) on h and flu were used in the fitting procedure. 
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Fig. 4. — Root mean square deviation of the probability distribution functions of unlensed, 
lensed, and corrected standard sirens as a function of redshift. 
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Table 1. A summary of the results of the cosmology fits shown in Figure [2j Quoted 68.3% 
confidence level intervals for flu, wo, and h were obtained by finding the extremal value of 
the cosmological parameter 9 where log L{6) = logL max — 2.30/2. 



Simulated Data Set 


!2m 


Wo 


/( 


High-z 


Unlensed 


[0.25,0.32] 


[-1.6, 


-0.7] 


[0.64,0.82] 




Lensed 


< 0.40 


< - 


0.2 


> 0.57 




Corrected 


< 0.34 


< - 


0.4 


> 0.61 


High-z + low-z 


Unlensed 


[0.28,0.31] 


[-1.2, 


-0.8] 


[0.67.0.75] 




Lensed 


[0.20,0.38] 


[-1.6, 


-0.5] 


[0.63.0.80] 




Corrected 


[0.26,0.34] 


[-1.5, 


-0.6] 


[0.64.0.78] 


High-z + h known 


Unlensed 


[0.29,0.31] 


[-1.02, 


-0.98] 






Lensed 


[0.25,0.35] 


[-1.27, 


-0.83] 






Corrected 


[0.28,0.33] 


[-1.14, 


-0.94] 




High-z + Qm known 


Unlensed 




[-1.12, 


-0.79] 


[0.67,0.72] 




Lensed 




[-1.67, 


-0.38] 


[0.59,0.77] 




Corrected 




[-1.34, 


-0.56] 


[0.63,0.74] 


High-z + wo known 


Unlensed 


[0.29,0.31] 






[0.696,0.704] 




Lensed 


[0.25,0.37] 






[0.66,0.73] 




Corrected 


[0.28, 0.34] 






[0.68,0.71] 
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Table 2. A summary of the results of the cosmology fits shown in Figure [3j Quoted 68.3% 
confidence level intervals for wo and w a were obtained by finding the extremal value of the 
cosmological parameter 9 where log L{6) = logL max — 2.30/2. 



Simulated Data Set wq w& 



High-z 


Unlensed 


[-1.13, 


-0.85] 


[-0.6,0.4] 




Lensed 


[-1.41, 


-0.44] 


[-3.8, 1.5] 




Corrected 


[-1.29, 


-0.83] 


[-1.6,0.9] 


High-z - 


- low-z Unlensed 


[-1.13, 


-0.86] 


[-0.6,0.4] 




Lensed 


[-1.40, 


-0.59] 


[-2.8, 1.5] 




Corrected 


[-1.25, 


-0.83] 


[-1.4,0.9] 


High-z - 


- h and !)m known Unlensed 


[-1.04, 


-0.96] 


[-0.2,0.2] 




Lensed 


[-1.36, 


-0.69] 


[-1.6, 1.6] 




Corrected 


[-1.19, 


-0.86] 


[-0.8,0.8] 



Note. — In the fitting procedure Gaussian priors (3% accuracy) on h 
and Cm were used. 



